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Abstract 

We analyze supergravity models that predict a low mass gluino within the landscape of sparticle mass 
hierarchies. The analysis includes a broad class of models that arise in minimal and in nonminimal 
supergravity unified frameworks and in extended models with additional UiVfx hidden sector gauge 
symmetries. Gluino masses in the range (350 — 700) GeV are investigated. Masses in this range are 
promising for early discovery at the LHC at y/s — 7 TeV (LHC-7) . The models exhibit a wide dispersion in 
the gaugino-Higgsino eigencontent of their LSPs and in their associated sparticle mass spectra. A signature 
analysis is carried out and the prominent discovery channels for the models are identified with most models 
needing only ~ 1 fb _1 for discovery at LHC-7. In addition, significant variations in the discovery capability 
of the low mass gluino models are observed for models in which the gluino masses are of comparable 
size due to the mass splittings in different models and the relative position of the light gluino within the 
various sparticle mass hierarchies. The models are consistent with the current stringent bounds from the 
Fermi-LAT, CDMS-II, XENON100, and EDELWEISS-2 experiments. A subclass of these models, which 
include a mixed-wino LSP and a Higgsino LSP, are also shown to accommodate the positron excess seen in 
the PAMELA satellite experiment. 
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I. INTRODUCTION 



The exploration of supersymmetry at colliders within the landscape of sparticle mass hierarchies 
can provide insight into the nature of the underlying theory[l-^] . Thus, in the minimal supersym- 
metric standard model there are 32 supersymmetric particles including the Higgs bosons and in 
general there exists a mass hierarchy among them. One can estimate the number of sparticle mass 



patterns that may appear by use of Stirling's formula (n! ~ \/2irn(n/e) n ), which gives close to 10 28 
different possibilities including sum rules for the sparticle masses. Limiting oneself to the first four 
lightest sparticles and assuming that the lightest supersymmetric particle (LSP) is the neutralino, 
the number is significantly reduced but is still around O(10 4 ). On the other hand, in supergravity 
unified models with radiative breaking of electroweak symmetry, the number of allowed possibilities 
is much smaller [l-^]]. Further, it has been demonstrated in several subsequent works [2-8], that 
analyses of supersymmetry from the point of view of the possible low mass sparticle hierarchies 
indeed shed light on the underlying theory as the specific mass orderings and mass gaps among 
sparticles dictate the type of signatures that would be produced at the Large Hadron Collider. 
Specifically, in this work we investigate the landscape under the constraint that the gluino has a 
low mass and is discoverable in early runs at the LHC at ^Js = 7 TeV and with a few fb -1 of 

Q 

integrated luminosity (for a recent review see [9|). 

Because of their strong color interactions, gluinos and squarks are expected to be produced 
copiously in high energy pp collisions at the LHC, and a wide class of supergravity models are 
capable of such sparticle production. Several works have already appeared which point to encour- 
aging results for the possibility of discovery of supersymmetry in the early data 10h17I]. Here, 



our analysis of sparticle mass hierarchies is done within the context of supergravity grand unified 



models 



18[ where we include both universal and nonuniversal soft breaking 



20J] that 



allow for the possibility of low mass gluinos in the mass range of approximately (350 — 700) GeV, 
with sfermions that could be either light or heavy, consistent with all current experimental con- 
straints. In addition, we study extended supergr avity models where the relic density is satisfied by 



coannihilations with matter in the hidden sector 
solution, beyond a Breit-Wigner Enhancement 



23 



21 



22] which is another possible particle physics 



241 ] and other mechanisms, to achieving a large 



flux of dark matter in the halo consistent with the relic density. Dark matter annihilations in a 
universe with a nonthermal cosmological history is a possible solution as well 25|, 



The class of models we consider within the context of sparticle mass hierarchies are those which 
are promising for early SUSY discovery. A number of models with low mass gluinos and their 
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hierarchical mass patterns are exhibited, and their LHC signatures at yfs = 7 TeV are analyzed 
and discussed. The strong correlations between the hierarchical structure of the sparticle masses 
[l| and the early discovery prospects at LHC as well as in dark matter experiments are studied 
with specific focus on the gluino in the mass hierarchy. It is shown that most of the models would 
become visible with (1 — 5) fb _1 of integrated luminosity. Further, the models are subjected to 



the more recent constraints from the Fermi Large Area Telescope (Fermi-LAT) data 27J on the 
monochromatic radiation that can arise in annihilation of neutralinos into photons, as well as from 
the direct detection exp eriments on the LSP-nucleon spin-independent cross sections 
recent overview see 



28-laO 



30] (for a 



3l|) which are beginning to put more stringent constraints on super symmetric 
models. It is also shown that a subclass of the models with a low mass gluino have a neutralino 
LSP that can satisfy the positron excess as seen in the PAMELA satellite experiment 32| . The 



low mass gluino (LG) models discussed are consistent with the cold dark matter relic density from 
the WMAP data |33], and most lie within the observed WMAP experimental band. The sparticle 
mass hierarchies corresponding to these LG models are given in Table f|TT| ) . 

The outline of the rest of the paper is as follows: In Sec. (II) we discuss the experimental con- 
straints on the model classes studied. These constraints include the recent limits from Fermi-LAT 
on the monochromatic radiation and the upper limit on the spin-independent LSP-nucleon elastic 
scattering cross section. In Sec. (Ill) we analyze a number of representative model points that 
encompass a wide range of theoretical frameworks which include minimal supergravity grand uni- 
fication (mSUGRA), nonuniversal supergravity grand unified models (NUSUGRA), and extended 
supergravity models with a U{1)\ hidden sector gauge symmetry. The specific properties of these 
models, including their light sparticle spectra, are also discussed in this section. A subclass of the 
models is shown to be capable of accommodating the PAMELA data. In Sec. (IV) we discuss the 
signatures of the models considered here at the LHC with -y/s = 7 TeV. Conclusions are given in 
Sec.(V). 

II. EXPERIMENTAL CONSTRAINTS 

Below we summarize current constraints on supersymmetric models from collider and dark 
matter experiments which we include in the analysis. 

1. Collider Constraints: The models we consider are subject to several accelerator constraints 
which are sensitive to new physics. These include the experimental result Br{B — > X s ~f) = (352 ± 



23 ±9) x 10 6 from the Heavy Flavor Averaging Group 



3J along with BABAR, Belle, and CLEO, 
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while for the standard model prediction we use the next to next to leading order (NNLO) estimate 
of Br(b — > S7) = (3.15 ± 0.23) x 10 -4 [37]]. Supersymmetry makes important corrections to this 



process 



381 ]. and the current (small) discrepancy between the experiment and the standard model 



predictions hints at the presence of low- lying sparticles [40j] . For this analysis we take a 3a corridor 
around the experimental value: 2.77 x 10 -4 < Br(b — > S7) < 4.27 x 10 -4 . Another important 
flavor constraint is the rare decay process B s — > fJ- + fJ-~ ■ We take the recent 95% C.L. constraint by 
CDF [39] Br(B s < 5.8 x 1CT 8 . In the context of a large spin-independent cross section in 

the minimal super symmetric standard model (MSSM), which can occur at low Higgs masses, and 
large tan /3, this constraint has recently been shown to be very restrictive on the MSSM parameter 
space under radiative breaking Q], including models with low mass dark matter 4l| in the 10 GeV 
range The correction to — 2 is an important indication of new physics, and in models of 



44J. The 



supersymmetry, significant corrections are expected with a low-lying sparticle spectrum 

2 data [45J has been analyzed recently using improved estimates of the hadronic correction 
In this analysis we take a conservative bound: —11.4 x 10 -10 < 5(g^ — 2) < 9.4 x 10~ 9 . 
Experimental limits on the sparticle masses are as discussed in {40! ]. 

2. Dark Matter Constraints: In addition to the above, there are constraints from dark matter 
experiments. The seven-year WMAP data [33|] gives the relic density of cold dark matter as 
Qdm'i 2 = 0.1109±0.0056. For our analysis, we take a conservative upper bound on the relic density 
predictions in the MSSM to take into account the theoretical uncertainties and other possible dark 
matter contributions (see e.g., j^j], [43]], [48]; for early work see 49]). However many of the models 
produce a relic abundance consistent with the double-sided WMAP bound at a few standard 
deviations. Further, new experimental results from probes of direct and indirect detection of dark 
matter have begun to place more stringent bounds on dark matter models. Thus, recent results 
from direct detection experiments CDMS-II 28], XENON100 [29], and EDELWEISS-2 Q suggest 



an upper bound to the spin-independent (SI) LSP-nucleon scattering cross section. At present, the 
limit is ctsi < 5 x 10~ 8 pb over the range of interest of the models discussed here. Many of the 
models discussed here have spin-independent cross sections in the range < (10 -9 — 10 -8 ) pb and 
are thus directly relevant to dark matter direct detection experiments in the near future. 

Further, the Fermi-LAT experiment has obtained constraints on the production of 7-ray lines 
with energies from 30 GeV to 200 GeV. The upper limits of the 7-ray line flux are in the range of 
(0.6— 4.5) x 10 -9 cm _2 s _1 [27]]. This gives rise to upper bounds on the corresponding dark matter 
annihilation cross sections, and as we will discuss later, the data constrains the eigencontent of 
the LSP in the mass range analyzed. Finally, the recent PAMELA 32| data shows an anomalous 
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A sample of low mass gluino models within the landscape 



Label 


M- g 


mo 


M 1 


M 2 


M 3 


A 


tan/3 


LG1 


424 


2000 


130 


130 


130 


-1000 


8 


LG2 


715 


60 


300 


300 


300 


-100 


8 


LG3 


386 


1400 


800 


528 


132 


3000 


25 


LG4 


378 


3785 


836 


508 


98 


-6713 


20 


LG5 


385 


2223 


859 


843 


130 


4680 


48 


LG6 


391 


1180 


860 


790 


138 


2692 


42 


LG7 


442 


2919 


263 


151 


138 


4206 


18 


LG8 


417 


1303 


257 


152 


139 


1433 


18 


LG9 


696 


1845 


327 


193 


249 


1898 


13 


LG10 


365 


1500 


1600 


1080 


120 


2100 


15 


LG11 


433 


605 


302 


176 


161 


1781 


22 


LG12 


588 


636 


419 


249 


228 


1568 


37 


LG13 


684 


335 


391 


466 


279 


-1036 


3 


LG14 


618 


48 


289 


310 


256 


-407 


6 


LG15 


602 


61 


310 


351 


249 





9 


LG16 


343 


2200 


450 


235 


100 


300 


5 


LG17 


425 


3000 


400 


207 


125 
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TABLE I: A sample of low mass gluino models where additionally we take (i > and m t ( po [ e ) = 173.1 GeV. 
The soft breaking parameters are given at the high scale of ~ 2 x 10 16 GeV . Nonuniversalites in the gaugino 
sector M 0= i j2 ,3 are taken in 15 of the models. All masses in the table are in GeV. A more detailed discussion 
of the models is given in the text. 

high positron flux in the range (10— 100) GeV. The possible dominant contribution to the positron 
excess in the MSSM with neutralino dark matter can arise from XiXi ~^ W + W~ , ZZ annihilations, 
depending on the eigencontent of the LSP. Here, a velocity averaged cross section on the order of 
(crv)ww+zz ^ 5 x 10~ 25 cm 3 s _1 will be shown to provide a reasonable explanation of the data 
without invoking large astrophysical boost factors in the positron flux. We will consider the 
implications of such constraints later in the analysis. 

III. MODELS OF A LOW MASS GLUINO AND SPARTICLE MASS HIERARCHIES 

In N = 1 supergravity unified models, gaugino masses can arise from the gauge kinetic function 
f a corresponding to the gauge group G a . Thus, for the standard model gauge groups SU(3)c, 
SU(2)l, and U(l) the gaugino masses are given by M a = (21Z(f a ))~ 1 F I dif a (a = 1,2,3), where 
the gauge kinetic function f a depends on the fields and F 1 are the order parameters for su- 
persymmetry breaking with the VEV of F 1 proportional to the gravitino mass, i.e., oc m 3 / 2 Mp. 
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A sample of the sparticle landscape for low mass gluino SUGRA models 



Label 


Mass Pattern 


Mg (GeV) 


Lightest Mg (GeV) 


Gluino Position 


LG1 


X? < xt < X°2 < 9 < X°3 < xl 


424 


1985 


4 


LG2 


x° < n < i R < v r < t> t < I L 


715 


635 


31 


LG3 


Xi < 9 < xt < X°2 < ii < xl 


386 


1411 


2 


LG4 


Xi < 9 < xt < X° < h < Xs 


378 


3751 


2 


LG5 


Xi < 9 < xt < X° < h < fi 


385 


2217 


2 


LG6 


X? < h < g < xt < X°2 < Xs 


391 


1202 


3 


LG7 


Xi < xt < X°2 < 9 < X°3 < Xa 


442 


2888 


4 


LG8 


Xi < xt < X°2 < 9 < X°3 < Xl 


417 


1314 


4 


LG9 


Xi < Xi < xl < xi < xl < xt 


696 


1882 


7 


LG10 


Xi < Xi < X° < 9 < Xs < *i 


365 


1511 


4 


LG11 


Xi < Xi ^ X°2 < ti <9<li 


433 


690 


5 


LG12 


X? < xt < Xa < Tt < Xs < X4 


588 


790 


14 


LG13 


x? < k < xt < xl < h < Ir 


684 


677 


24 


LG14 


Xi < h < Ir < v T < VI < h 


618 


550 


31 


LG15 


Xi < fi < Ir < v T < v t < xt 


602 


536 


31 


LG16 


x? < xf < xl < 9 < xs < xi 


343 


2178 


4 


LG17 


Xi < xt < X2 < 9 < Xs < xt 


425 


2966 


4 



TABLE II: Sparticle mass hierarchies for the low mass gluino models. Listed are the first six lightest 
sparticles in the spectra. The lightest squark, shown in column 4, is taken from the first 2 generations. In 
the display of the mass hierarchies we do not include the light CP even Higgs. 

The assumption that the fields cj>i are singlets of the standard model gauge group will lead to 
universal gaugino masses. However, in general the breaking can occur involving both singlet and 
non-singlet terms (for a recent analysis see jsol ]) and, thus, the gaugino masses in general will not 
be universal. Additionally, there can be loop corrections to the gaugino masses via exchange of 
heavy fields where the leading corrections are also proportional to 7713/9, to the Casimir and to the 
square of the gauge coupling for each gauge group (for early work see [5 11] ). 

A case of immediate interest for LHC analyses is M^{Mu) £ (100 — 300)GeV as the gluino 
mass at the unified scale (Mjj) increases in magnitude as one moves towards the electroweak 
scale via the renormalization group equations, and this translates to gluino masses in the range of 
(350 — 750) GeV in the model classes considered. While the models considered in this analysis can 
be obtained with specific choices of supersymmetry breaking, involving, for example, a combination 
of singlet and non-singlet F-term breaking for gaugino masses at the scale Mjj, we will keep our 
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analysis rather generic in that the information on symmetry breaking is contained in the gaugino 
masses at the GUT scale which can arise from various models of the gauge kinetic energy function, 
and in general they depend on both the hidden and the visible sector fields. 

Indeed there are several classes of well motivated SUGRA models which lead to sparticle spectra 
with a low mass gluino. These include the mSUGRA models and nonuniversal SUGRA models 
where the soft breaking in the gaugino sector at the grand unification scale includes non-singlet 
contributions as well as singlet and non-singlet contributions (see, e.g., 



50 



520). Different high 



scale models can generate different next to lightest supersymmetric particles (NLSPs) at the elec- 
troweak scale, some of which lead to the gluino as the NLSP (GNLSP models) uncovered in Ref. 



and further investigated in 



3-. 



5214541] . We also discuss models where the relic density due to ther- 



mal annihilations is rather small but can be boosted to the current range of the WMAP value with 
coannihilations with matter in the hidden sector 1 221. Such a situation arises in SUGRA models 



with extended U(l)x gauge groups 22]. Another possibility to boost the relic density arises from 



nonthermal processes 
ries, namely, thermal 



17 



Both of these mode 



22 1 and nonthermal 



25. 



asses, with different types of cosmological histo- 



261 ] . can lead to an explanation of the positron 



excess seen by the PAMELA satellite experiment and provide explanations for the relic abundance 
of dark matter - however, both classes require fields beyond those in the low energy MSSM. 

We give a broad sample of models in Table(jl]), where the SUGRA model parameters at the GUT 
scale are exhibited. The parameters include the universal scalar mass mo, the three gaugino masses 
-^l,2 3> the universal trilinear coupling Aq, and the ratio of the two Higgs vacuum expectation 
values, tan f3. These models satisfy the stringent Fermi-LAT constraints on 77 and 7Z cross 
sections as discussed later. The sparticle mass hierarchies corresponding to these models are given 
in Table (jH"| ), In Fig. ([I]) we give a display of the significance for several optimal channels for a 
sample of the model classes at 1 fb _1 of integrated luminosity at yfs = 7 TeV. One finds that 
all of the models shown are visible in several channel except for LG3 and LG4, which are visible 
only in a small number of channels. We discuss these results in depth in Sec. (IV). A more detailed 
discussion of these model classes follows. 

1. Low mass gluinos in mSUGRA: The possibility that a low mass gluino and heavy scalars 
can arise in the radiative breaking of the electroweak symmetry (REWSB) in SUGRA models was 
seen early on [55j]. It was later realized that this phenomenon is more general, and the Hyperbolic 
Branch (HB) or the Focus Point (FP) region of REWSB was discovered where scalars are heavy 



and gauginos are light 



34j. In the analysis of Q], it was found that on the HB the chargino is 



7 



50 



40- 



g 30 

CO 

o 
"E 

.5> 20 
CO 



10- 









LHC 7 TeV 



L = 1 fb 



-1 





LG1 LG3 LG4 



LG8 LG10 LG11 
Models 



LG14 LG17 



FIG. 1: (Color online) A summary of the significance (SUSY signal divided by the square root of the SM 
background) for various signature channels/cuts (C) for a subset of the low gluino (LG) mass models with an 
integrated luminosity of 1 fb -1 at the LHC for y/s — 7 TeV. We will discuss in detail in Sec. (IV) how these 
drastic variations come about even though the gluino masses here are confined to the range ~ (350-700) 
GeV for all models. 



predominantly the NLSP over a very broad class of soft breaking models . 

However, a low mass gluino, around 400 GeV, can arise in minimal supergravity if the neutralinos 
annihilate near the Z-pole or near the Higgs poles 



561 ], Specifically a low mass gluino arises in 
mSUGRA in the sparticle mass landscape which was labeled as mSP4 in [1] where the first four 
sparticles have the mass hierarchy as follows: 



Model LG1: X i < Xx < X° < 9, M~ g = 424 GeV. (1) 



An illustration of such a model is given by LG1 in Table (|IIip . In this case the LSP is essentially 
dominated by the bino component 2 as can be seen from Table (IIIip . The model LG1 has a low 



1 HB models with larger Higgsino components and with a light gluino have been studied in [34j-t3t|,[lb-3, 52], and 
Ref. 1 of 0- 

2 We use the notation Zn, Z\2, Z13, Z14 for the eigen components of the lightest neutralino, namely the bino, the 



wino, and the two Higgsinos. 
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Low mass gluino models in minimal SUGRA 



Label 


M h 






M S 


M fl 




M A 




Zu 


Z\2 


Zl3 


Z14 


er?o (cm 2 ) 


LG1 


115 


53 


104 


424 


1983 


1118 


2040 


1610 


0.993 


-0.051 


0.103 


-0.025 


5.4 x 10~ 46 


LG2 


111 


118 


220 


715 


126 


476 


463 


606 


0.990 


-0.044 


0.127 


-0.051 


5.2 x 10~ 45 



TABLE III: mSUGRA models of a low mass gluino and a display of some of the lighter masses within the 
sparticle mass hierarchies. In model LG1 the neutralino, the chargino and the gluino are all light while the 
SUSY scalars are heavy. This model generates the relic density in the WMAP band via the annihilation of 
the neutralinos near the Higgs pole. However, in model LG2 the scalars are also light and the relic density 
lies in the WMAP band via coannihilations of the neutralinos with the stau and other light sleptons. All 
masses are in GcV. 

mass gluino of about 420 GeV, and the relic density is within the WMAP band due to neutralino 

annihilations near the Higgs pole. It is easily seen that the model LG1 is on the HB and, thus, has 
reavy scalars. Since the neutralino is bino-like this leads to scaling of the gaugino mass spectra 
5t| . Further, because the neutralino has a suppressed Higgsino content, the spin-independent 

scattering is suppressed at the level of ~ 5 x 10 -46 cm 2 . 

A diametrical situation holds for the model point LG2. While the gluino is still light in this 

model, lying in the sub-TeV region, it is actually the heaviest sparticle in the whole sparticle 

spectrum. Specifically the hierarchy is 



Model LG2 : x? < h < £r < v T < . . . < g, M~ g = 715 GeV. 



(2) 



The relic density in the model LG2 is satisfied via coannihilations with the stau and with 
the other sleptons. Since model LG2 has a lighter scalar sparticle spectrum, it has a larger spin- 
independent cross section than LG1 by a factor of about 10. The ordering of the gluino mass in 
the hierarchy will have dramatic effects on the signals at the LHC which we will discuss in detail 
shortly. 

2. Low Mass Gluinos in Nonuniversal SUGRA models: Low mass gluinos can arise from models 
where the gaugino masses are in general nonuniversal at the GUT scale, as for instance, from a 
combination of singlet and non-singlet F term breaking as discussed in the beginning of this section. 
Models LG3-LG17 in Table([T|) fall in this class. Such models can lead to a low mass gluino that can 
be the NLSP, the NNLSP, etc. The case when the gluino is the NLSP requires special attention 
since here the relic density can be satisfied by the neutralino coannihilations with the gluino, as 
in the model LG3. This case was discussed at length in 52J, 



54, 3- 



For the models discussed 



in Table (|IVp . which have a gluino NLSP (GNLSP) with a rather low mass gluino (i.e., below 400 
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Low mass gluinos in GNLSP models 



Label 


M h 


M xl 


M xt 


M S 


M fl 


ti 


M A 




Zn 


Z\2 


Zia 


Z14 


cr?o (cm 2 ) 

XiP v ; 


LG3 


112 


340 


429 


386 


1253 


455 


1421 


995 


0.997 


-0.026 


0.067 


-0.029 


6.3 x 10~ 46 


LG4 


125 


377 


454 


378 


3529 


1244 


3888 


2615 


0.999 


-0.006 


0.021 


-0.005 


8.4 x 10~ 48 


LG5 


117 


365 


660 


385 


1081 


679 


1167 


1321 


0.999 


-0.003 


0.039 


-0.012 


2.1 x 10~ 46 



TABLE IV: A display of the lighter sparticle masses within the mass hierarchies, and other attributes of 
GNLSP models with low mass gluinos. The mass splitting between the gluino and neutralino is between 
~ (1 — 50)GeV for these models. Further details are given in the text. All masses are in GcV. 



GeV), the sparticle mass hierarchies are given by 52| 



Models (LG3, LG4, LG5) : Xi < 9 < Xf < X2 < h, M~ g = (386, 378, 385) GeV 



(3) 



Since the gluino coannihilation in the GNLSP model is a relatively new entry among the ways 
dark matter originates in the early Universe, we summarize the main features of the relic density 
calculation here first, before discussing the relevant features of this class of models that enter in the 
LHC signature analysis. Thus, for a GNLSP the relic density depends strongly on coannihilation 
effects which are controlled by the Boltzmann factor 58] 

9i (l + A^e-^ 



eq 

n- 



(4) 



™ eq J2 j 9j( 1 + A j) 3/2e ~ AjX ' 
where gi are the degrees of freedom of Xii x = m i/T, and Aj = (rrij — mi) /mi, with mi defined 
as the LSP mass. Thus, for the analysis of the relic density, the effective annihilation cross section 
cr e ff can be written approximately as 



(5) 



where we have used the fact that the gluino annihilation cross sections are usually much larger 
than the LSP annihilation even with inclusion of the Boltzmann factor and 

7TCC/3 



<r(gg^qq) = ^^(3 - /3 2 )(3 - 
lops 

l~~ \ c 37ra s 



log i±| [21 - 6/3 2 - 3/3 4 ] - 33/3 + 17/3 3 } , 



(6) 
(7) 



where the non-perturbative corrections to the annihilation cross section can arise via multiple gluon 
exchange, giving rise to a Sommerfeld enhancement factor £ . These effects may be approximated 



by 



60] 



Cjira s 



1 — exp 



CjTTCts 



(8) 
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where Cj =g = 1/2(C, =9 = 3/2) for gg — > gg(gg — > qq). In the above (3 = wl — 4m|/ 's, and 



/3 = ^1 ~~ 4m^/s. Although the gluino annihilation cross section a gg varies with gluino mass, the 
Boltzmann suppression factor 7^ controls the contribution to the <r e g so that the relic density of 
the bino-like neutralino is consistent with WMAP. We note that the effects of the Sommerfeld 
enhancement on the gluino cross sections can increase by a small amount for a bino-like LSP 



~ (2 — 3)% [52J, [59J]. Such an increase in the mass gap between the gluino and the neutralino can 
potentially enhance the discovery reach of this class of models at the LHC as the mass gap between 
the g and Xi pl avs a crucial role in the strength of the LHC signals [52] which we will discuss in 
details later. Three GNLSP models are exhibited in Table (|IVp . Some of their pertinent spectra 
and other attributes, including their spin-independent cross sections, are also given in Table (|IV|) . 
It is seen from this table the neutralino is dominantly a bino. For each of these models gluino 
coannihilation dominates the relic density calculations. In the absence of additional hidden sector 
gauge groups (to be discussed in what follows), only LG3 lies in the WMAP band, while models 
LG4 and LG5 each have a reduced relic abundance and reduced mass gaps between g and x?- 

Next we consider five models, LG6-LG10, in Table([T]) where the gluino is light and in some 
cases it is the next to next to LSP (GNNLSP). Specifically the sparticle mass hierarchies are 

Model LG6 : X°i < h < 9 < xf < X% 

Model LG7 : Xi < Xi < X2 < 9 < xl 

Model LG8 : X°i < xf < X° 2 < 9 < xl 

Model LG9 : x? < xf < xl < Xs.4 < xt < 

Model LG10: x? < xf < X% < 9 < xl < *l, M g = 365 GeV . (9) 

Their relevant sparticle spectra and other attributes including eigencontent, and spin-independent 
cross sections are displayed in Table (fV] ). Here the NLSP can be the stop or the chargino, while the 
LSP is a mixture of bino, wino and Higgsino. The specific nature of the neutralino eigencontent 
makes these models significantly different from each other and from other light gluino models. This 
includes an interesting subclass of models where the LSP is dominantly Higgsino (LG10). Such a 
model class was recently analyzed in the context of both the PAMELA positron excess and 
the Fermi-LAT photon line data. 

Finally, in Table (]Vip we give five models, LG11-LG15 where in addition to the low mass gluino 
one also has a light stau and a light stop (due to the smaller GUT value of mo) along with a light 
chargino. Models LG11-LG15 have compressed sparticle spectra with the heaviest sparticle mass 



M 9 


= 391 


GeV , 


M~ g 


= 442 


GeV , 


M g 


= 417 


GeV , 


M g 


= 696 


GeV , 


M g 


= 365 


GeV . 
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Low mass gluinos including GNNLSP models 



Label 


M h 


M x1 


M xt 


M g 






M A 




Zu 


Zl2 


Z13 


Zl4 


cr?o (cm 2 ) 

XiP K ' 


LG6 


109 


359 


570 


391 


737 


378 


839 


833 


0.992 


-0.020 


0.109 


-0.064 


6.6 x 10~ 45 


LG7 


119 


108 


120 


442 


2786 


1357 


2947 


2185 


0.998 


-0.049 


0.039 


-0.006 


2.4 x 10~ 47 


LG8 


112 


104 


117 


417 


1255 


718 


1284 


1032 


0.968 


-0.202 


0.143 


-0.041 


2.7 x 10~ 45 


LG9 


115 


135 


152 


696 


1811 


1067 


1874 


1513 


0.985 


-0.136 


0.098 


-0.030 


7.8 x 10~ 46 


LG10 


112 


111 


115 


365 


1570 


734 


1609 


1323 


0.058 


-0.075 


0.721 


-0.686 


7.3 x 10~ 45 



TABLE V: The spectrum of low mass sparticles including the GNNLSP models within the sparticle mass 
hierarchies, and other attributes of models in NUSUGRA with a low mass gluino. Model LG6 is a GNNLSP, 
and models LG7,LG8, and LG9 are effectively GNNLSP as the chargino and second heaviest neutralino are 
roughly mass degenerate. LG10 has a mass splitting between the chargino and the neutralino of <~ 5 GeV 
and is effectively a GNNLSP model. All masses are in GeV. 



Models with a low mass gluino with a light stau and a light stop 



Label 


M h 


M x1 


M-± 

Xi 


M g 




ti 


M A 


M- bl 


Zu 


Z\2 


Z\z 


Z14 


cr?o (cm 2 ) 

XiP K ' 


LG11 


103 


121 


134 


433 


516 


271 


701 


489 


0.978 


-0.171 


0.117 


-0.032 


6.6 x 10~ 45 


LG12 


107 


169 


190 


588 


424 


489 


551 


602 


0.973 


-0.171 


0.146 


-0.058 


3.0 x 10" 44 


LG13 


104 


160 


359 


684 


365 


187 


845 


604 


0.996 


-0.023 


0.076 


-0.041 


2.0 x 10~ 45 


LG14 


111 


114 


227 


618 


118 


338 


478 


514 


0.990 


-0.043 


0.124 


-0.052 


4.7 x 10~ 45 


LG15 


109 


121 


238 


602 


132 


397 


398 


525 


0.979 


-0.052 


0.177 


-0.083 


1.7 x 10" 44 



TABLE VI: An exhibition of the light sparticle masses within the hierarchies in models with low mass 
gluinos that are accompanied by light staus and light stops. The light Higgs masses for these models lie in 
the range (103-111) GeV and are generally lighter than the models where the chargino or the gluino is the 
NLSP shown in the previous tables, as the scalars in such model are much heavier. All masses are in GeV. 



around 850 GeV. Model LG11 has a highly reduced overall mass scale of the sparticles, with a low 
mass gluino and a light stop with the mass hierarchy 

Model LG11 : x°i < xt < X°2 < h < g, M g = 433 GeV , (10) 

with the remaining sparticles in the mass range (500-700) GeV. In each of these models the relic 
density can be satisfied via coannihilations with different superparticles, and, in particular, model 
LG13 proceeds via stop coannihilations. For LG12, the gluino mass lies in the middle of the 
sparticle mass spectra, while for (LG13-LG15) the gluino mass is close to being the largest mass 
even though it is still relatively light, i.e., M g < 700 GeV. Models LG11 and LG13 have rather low- 
lying light Higgs. However, the extraction of the Higgs mass from LEP data is model dependent, 
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Low mass gluino SUGRA models which explain the PAMELA positron excess 



Label 


M h 




M xt 


M S 




ti 


M A 


M- hl 


Zn 


Z\2 


Zl3 


Zl4 


cr?o (cm 2 ) 

XiP v ; 


LG10 


112 


111 


115 


365 


1570 


734 


1609 


1323 


0.058 


-0.075 


0.721 


-0.686 


7.3 x 10~ 45 


LG16 


112 


182 


188 


343 


2193 


1254 


2290 


1780 


0.654 


-0.718 


0.209 


-0.112 


3.0 x 10~ 44 


LG17 


111 


168 


171 


425 


2986 


1698 


3215 


2421 


0.724 


-0.681 


0.101 


-0.043 


5.9 x 10~ 45 



TABLE VII: An exhibition of the sparticle mass hierarchies in low mass gluino models which can explain 
the PAMELA positron excess. The models are consistent with data from the Fermi-LAT [23, CDMS-II 
01, XENON100 and EDELWEISS-2 Q experiments. In the above we do not list the component of 
the neutralino lying in the hidden sector as it is typically small, i.e., Z\ k < 1%. 



and we retain these models in the analysis pending further experimental data. 

3. Low mass gluinos in extended SUGRA models and PAMELA data: Recently, the PAMELA 
collaboration [32| has included an analysis of statistical uncertainties in the positron fraction 61] 



and presented its results on the absolute p flux [62]. Models LG10, LG16, and LG17 in Table(jT]) have 
low mass gluinos with many other desirable features. Specifically, they can explain the positron 



excess in the PAMELA satellite experiment [321 ] (for previous experiments see 63]). However, 
their relic density by the usual thermal annihilation processes is rather small. This situation can 
be modified at least in two ways. 

First, the relic density could be enhanced via coannihilation with matter in the hidden sector. 
For example, an extended Abelian gauge symmetry can arise from the hidden sector and can 
couple to the MSSM sector via mass mixing or kinetic mixing with the hypercharge field. The 
mass mixing arises via the Stueckelberg mechanism and the kinetic mixing via mixing of the 
hidden sector field strength and the hypercharge field strength. Such mixings can lead to an 
enhancement of the relic density which we now discuss. Thus, we consider for specificity the 
mass growth via the Stueckelberg mechanism allowing for mixings between the hypercharge gauge 
multiplet (Y^, Ay , Ay , Dy) and the gauge multiplet of U{l)x (AT«, Ay, Ay , Dx) both taken in the 



Wess-Zumino gauge. The effective Lagrangian contains the mixing terms 
and Ref. 



2H 



64J, Ref (1-4) of [65(], 



-(d^a + MyY^ + MxX^f + [Vst(M X A X + My Ay) + h.c] , 



(11) 



where the axionic field a arises from a chiral multiplet S = (p + io~, ip s t, F$). The fields ^ s t and Ax 
produce two Majorana spinors (hidden sector neutralinos) which mix with the neutralinos in the 
visible sector via the mass mixing given by Eq. (jll|) . In addition the hidden sector neutralinos can 
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mix with the visible sector via kinetic mixing 2l|, [66j 



- -X^Y^ - i5{\ x a ■ d\ Y + {Y o X)) + 5D X D Y . (12) 
We note that My : Mx and (5 are constrained to be small by fits to the precision electroweak data 



651 ] and thus the additional neutralino states are weakly coupled to the MSSM. 

An estimate of the bound on the (mass or kinetic) mixing can be obtained from studying 
eigenvalues of the vector sector. One finds the mixing is constrained for masses near the Z-pole as 
derived in Ref. 1 of 65] 



|e| <. 05^ 1 - ml/ml , (13) 

where e is the overall mass and/or kinetic mixing and J7i(i 2) is the mass of the hidden sector vector 
boson, or the Z boson mass, depending on which side of the Z-pole the hidden sector mass resides. 
We add that dark matter with Stueckelberg mass growth and/or kinetic mixing, with a massive 



hidden U(l) [211 ] has also appeared in Refs. 67|] in various contexts. 

The analysis above can be extended to a U(l) x gauge group in the hidden sector that produces 
2n additional Majorana fields in the hidden sector, which we denote collectively by A , (k = 1 — 2n). 
This extended set of hidden sector Major anas will mix with the visible sector via mass mixing and 
kinetic mixing, where we again take the mixings to be small. Due to the small mixings between the 
visible and the hidden sector Majoranas, the LSP in the visible sector will have a small component 
which lies in the hidden sector, and as such, the eigencontent of the LSP is modified so that 
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22, 



64 



68] 



2n 

X? = Z n B + Z 12 W + Z^Hr + Z U H 2 + £ , (14) 

k=l 

where are the hidden sector neutralinos (Stinos) composed of the hidden sector gauginos and 
the hidden sector chiral fields as discussed above and Z^ k record the leakage into the hidden 
sector. Because of the small mixings between the visible and the hidden sector, the Z^ k are rather 
small, typically less than 1% of the components in the visible sector. For this reason we do not 
record them in Table (|VIIj) and elsewhere. However, we will assume that the overall mixing is large 
enough so that the extra- weakly interacting states [^.[^ remain in contact with the thermal bath. 
Specifically, we envision mixings in the range (10~ 5 — 10 -2 ). In what follows, in the context of the 
PAMELA data, we will assume that the hidden sector states lie above the Z-pole. The system of 
additional Majoranas then mostly have hidden sector eigencontent and have masses close to the 
mass of the LSP. 
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NFW min, p = 0.40 GeV/cm 3 , t E = 2 

SUGRA models 



PAMELA Data 
PAMELA Data (stat) 




Mixed-Wino LG17 

Mixed-Wino LG16 

— Galprop Background Fit 



E GeV 




Kinetic Energy T GeV 



FIG. 2: (Color online) Top: Mixed-wino LSP models LG16 and LG17, the Higgsino LSP model LG10, 
and the PAMELA positron excess. For comparison a pure wino LSP model is also shown. Boost factors 
of (1-3) are used. We note that an LSP with a pure wino content while compatible with the PAMELA 
data has difficulties with the 77 and 7Z cross sections from Fermi monochromatic photon data. For this 
reason Tabled does not contain a pure wino model. Bottom: The PAMELA anti-proton flux and the SUSY 
models. The p flux is compatible with the data. The mixed-wino models as well as the Higgsino model can 
produce large LHC jet signatures as discussed in the text. 

We consider now the inclusion of the hidden sector which may have a very significant effect on 
the relic density. For the case that the Major anas in the hidden sector are effectively degenerate 
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in mass with the LSP, the LSP can coannihilate with the hidden sector neutralinos generating an 
enhancement for the density of relic neutralinos relative to that in the MSSM, through the extra 
degrees of freedom supplied by the hidden sector, so that 

n^oh 2 ~ f E x nff SM h 2 , (15) 



where Je is the relic density enhancement factor given by 

Je 
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^ dhid 1 " 



dL 



(16) 



and where cihid is the degree of degeneracy in the hidden sector and d v is is the degree of degeneracy 



in the visible sector. We note in passing that an analysis on t 



re effects of fields in the hidden sector 



on the relic density via coannihilation effects was studied in 2l| and a similar idea was pursued in 
[70I but without an additional hidden sector. For the J7(l)y extension of the SUGRA model with 
the hidden sector neutralinos essentially degenerate in mass with the LSP, one has cfhid = 2 x 2n, 
where 2n is the number of Majorana fields in the hidden sector with 2 degrees of freedom for each 
Majorana field. For the case when the LSP neutralino is pure wino, the chargino is essentially 
degenerate in mass with the LSP and in this case d v j s = 2(LSP) + 4(chargino), while for the 
case when the LSP neutralino is a mixture of bino, wino, and Higgsinos, one has d v ; s — > 2(LSP). 
Thus, for the case when the neutralino and the chargino are not degenerate as is the case when the 
eigencontent of the neutralino is a mixture of bino, wino and Higgsino, a large Je can be generated, 
for example, Je = 49 for n = 3, where the maximal Je we obtain is ~ 40 due to coannihilations 
in the visible sector. We will utilize this enhancement of the relic abundance in fitting the WMAP 
data. 

With the above enhancement mechanism one can achieve a large (<7t>)haio needed for a solution to 
the PAMELA data and at the same time have a relic density close to the WMAP data. This indeed 
is the case for models LG16 and LG17 which can explain the PAMELA data via the neutralino 
annihilation X1X1 ~^ W + W~ with the subsequent dominant decay W + — > £ + V£ while at the same 
time satisfying the relic density constraint via the additional hidden sector degrees of freedom. In 
these cases, the mixed-wino annihilation cross section into WW ~ 6x 10 -25 cm 3 s _1 , and only small 
effects from the possible clumps in the local halo are needed. For LG10, the Higgsino-like LSPs 
annihilate into both WW and ZZ, where in addition the decay Z — > £ + £~ gives a contribution to 
the positron excess. Here one finds ((tv)ww+zz ~ 3.5 x 10~ 25 cm 3 s~ 1 and one needs a clump factor 
of around (2-3) to fit the PAMELA data while the LSP mass is significantly lower, ~ 111 GeV 
[r?| . than the mixed wino case. Interestingly, this Higgsino-like model produces smaller p flux at 
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/ \ Theory j / \ Theory 

(av) lZ and (av)^ 



Label 


E 1 (GeV) 


( CrW )7Z C ° Iy ( Cm3 / S ) 


E 1 (GeV) 


(crw)™ oory (cm 3 /s) 


LG1 


14 


2.7 x lO" 30 


53 


1.5 x 10~ 34 


LG2 


101 


2.6 x 10- 30 


118 


9.4 x 10~ 30 


LG3 


334 


1.7 x 10~ 32 


340 


2.5 x 10~ 32 


LG4 


372 


2.2 x 10~ 34 


377 


6.9 x 10~ 34 


LG5 


359 


1.1 x 10~ 33 


365 


9.1 x 10~ 33 


LG6 


353 


3.2 x 10~ 32 


359 


6.7 x 10~ 32 


LG7 


89 


1.2 x 10" 28 


108 


2.8 x 10~ 29 


LG8 


84 


1.5 x 10~ 28 


104 


3.5 x 10~ 29 


LG9 


119 


8.8 x 10~ 29 


135 


1.7 x 10~ 29 


LG10 


92 


1.8 x 10" 28 


111 


8.7 x 10~ 29 


LG11 


103 


4.0 x 10" 29 


121 


7.1 x 10~ 30 


LG12 


157 


3.5 x 10~ 29 


169 


5.0 x 10~ 30 


LG13 


147 


1.3 x 10~ 31 


160 


4.0 x 10~ 31 


LG14 


96 


3.1 x 10~ 30 


114 


1.1 x 10~ 29 


LG15 


104 


3.1 x 10~ 30 


121 


7.9 x 10~ 30 


LG16 


171 


6.0 x 10~ 27 


182 


1.1 x 10~ 27 


LG17 


156 


9.4 x 10~ 27 


168 


1.7 x 10~ 27 



E 1 (GeV) 


NFW 


Einasto 


Isothermal 


Model MgSfi^ 


(90-100)[110-120] 
(170-180)[180-190] 
(150-160) [160-170] 


(6.0-3.8)[1.0-1.6] 
(4.0-6.1) [2.7-3.2] 
(8.2-6.3) [2.7-1.7] 


(4.3-2.8) [0.7-1.1] 
(2.9-4.4) [1.9-2.3] 
(5.9-4.5) [2.0-1.3] 


(10.3-6.6)[1. 7-2.7] 
(6.8-10.4) [4.6-5.5] 
(14.1-10.9) [4.7-3.0] 


LG10 (0.18)[0.087] 
LG16 (6.0) [1.1] 
LG17 (9.4)[1.7] 



TABLE VIII: Top table: Predictions for x?X? ~> 77 1 1% f° r the models including the bino-like models, 
Higgsino-likc models and wino-like models. Lower table: The current experimental constraints from Fcrmi- 
LAT [27| on 77 and 7Z modes are shown with three halo profiles. Constraints are shown for those model 
classes which can describe the PAMELA data, and have large neutralino self annihilation cross sections into 
77 and 7Z. Shown are the ranges where the models can be constrained where the notation is [..] for the 77 
mode and (..) is for the mode. All cross sections in the lower table are given in 10 _27 cm 3 s _1 . Models 
LG16, LG17 are within reach of the Fermi-LAT limits as they have a large wino content, and are discussed 
in detail in the text. 



higher energies and slightly larger flux at lower energies relative to the heavier mixed wino case. 
Also the Higgsino-like LSP receives less stringent constraints from the Fermi photon line data than 
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the wino-like LSP. We note that models with a Higgsino-like LSP have been studied in the past 
711 ] . However, to our knowledge the model class uncovered in [17j, for which LG10 is an example, 
is the first illustration of a SUGRA model with a Higgsino LSP that gives a fit to the PAMELA 
data. 

To elaborate further, the models of Table (|VIip have a neutralino mass in the range of (111 — 
183) GeV which allows a (crv)haio m the range (3.5 — 6) x 10~ 25 cm 3 s _1 , and thus the model 
can produce the positron excess seen in the PAMELA experiment. Further, as seen in Table !) VIip 
these models are consistent with the current XENON100 limits, and have the possibility to produce 
scattering cross sections that are discoverable in improved dark matter experiments, i.e., in the 
range ~ 5 x (10 -45 — 10~ 44 ) cm 2 . In addition, these models are consistent with the Fermi-LAT 
photon data as seen in Table tjVIIip . Hence, the eigencontent of the LSP in these models is mixed- 
wino for LG16 and LG17, and close to being pure Higgsino for LG10, and therefore their dark 
matter signatures are modified drastically relative to other models classes discussed which are 
significantly bino-like. 

In the top panel of Fig. ([2]) we show several fits to the PAMELA positron fraction where we 
include the data from recent experiments 



ment with 



32, 



611 ] . The analysis of this panel shows general agree- 



131 . 1 221 . |72[ ] for the pure or essentially pure wino case (for early work on cosmic rays 



relevant to this discussion see 



73j and [74I]). Here, however, we specifically show that the mixed- 



wino models LG16 and LG17 and the Higgsino LSP model case, LG10, provide a good description 
of the PAMELA data with boosts of order unity. For comparison we also show the essentially pure 
wino case, which requires no boost, but has difficulty explaining the Fermi photon data. In the 
lower panel of Fig.([2j) we give a comparison of the p flux with recently released data [(32 ]. Indeed it 
is seen that the theoretical prediction of the p flux is in good agreement with this data. The boost 
factors used are rather minimal and are not assumed different for the e and p fractions, which is in 
principle a possibility (different boosts are often introduced in order not to upset the p flux while 



enhancing the e flux). The analysis we present does not attempt to explain the AT 
energy e + e data. Such data could be ex plai ned with an additional electron source 
wino LSP 



2 1, 0, EJ or mixed-wino LSP 



13, 
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a 



^ermi high 
75| with a 



IV: SIGNATURE ANALYSIS AT THE LHC AT *Js = 7 TEV 



Signature analyses at center of mass energies of ^/s = 10 TeV and y/s = 14 TeV at the LHC 
already exist in the literature, and as mentioned in the introduction, a few analyses at the center 
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1CH171]. For this analysis, our emphasis is on 



of mass energies of y/s = 7 TeV have also appeared 
the discovery of models which admit low mass gluinos in early runs at the LHC consistent with 
dark matter interpretations for a neutralino LSP. 

1. Standard Model Background: The discovery of new physics requires an accurate deter- 
mination of the standard model (SM) background. The recent works of 
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121 ] have given an 



analysis of such backgrounds including 2 — > n processes at \fs = 7 TeV appropriate for pp col- 
lisions at the LHC. We use for our analysis the simulated SM background^ of [12| which was 



generated with MadGraph 4.4 [76J] for parton level processes, Pythia 6.4 [77( for hadronization 
and PGS-4 [T8f] for detector simulation. An MLM matching algorithm with a jet clustering 
scheme was used to prevent double counting of final states. Further, the 6-tagging efficiency in 
PGS-4 is based on the Technical Design Reports of ATLAS j3] (see 12]). which is similar to the 



efficiency of CMS 



801 ]. with the mis-tagging rate of 6-jet unmodified from the default in PGS-4. 



In addition Tauola is called for tau decays 



811 ] . The processes that are included in the SM back- 



ground are : (QCD 2,3,4 jets), (ii + 0, 1,2 jets), (66 + 0, 1,2 jets), (Z/7 (-> 11, vv) + 0, 1, 2, 3 jets), 
(W ± (-^ lv) + 0,1,2,3 jets), (Z/j(-> 11, vv) + it + 0, 1,2 jets), (Z/i(-> ll,vv) + 66 + 0, 1, 2 jets), 
(W ± (-»■ Iv) + 66 + 0, 1, 2 jets), (W ± (-> Iv) + it + 0, 1, 2 jets), (W ± (-> lv) + to (ib) + 0, 1, 2 jets), 
(ttti, ttbb, 6666), (W ± (-> lv) + W ± (-»• lv)), (W ± (-> lv) + Z (^ a//)), (Z (-)■ a//) + Z (-»• a//)), 
(7 + 1, 2, 3 jets). Here Z is e, t, a jet refers to gluon as well as first and second generation quarks 
and all denotes either Z, v or jet. The above processes are final state processes at the parton level, 
i.e. before hadronization. A more detailed discussion of the SM background can be found in 
which includes a list of cross section, number of events and luminosity for each process (see Table(I) 

of Q). 

2. SUSY Signal Generation and Optimization of signature cuts: The sparticle spectrum for the 
signal analysis was generated using SuSpect [82| via micrOMEGAs 



computed with SUSY-HIT 



831 ] and branching ratios are 



841 ]. Some differences in the output of the sparticle masses are known 



85] 



to exist when computed with different codes. We have checked our models with SOFTSUSY 
and found only small differences. 

We now discuss a set of signatures used in our early discovery analysis at the LHC at y/s = 
7 TeV. The notation used in the these signatures is: I denotes e,fi, and Pt(£), Pt{j) define the 
transverse momentum of the lepton t, and of the jet j, while n(£), n(j) give us the number of 
leptons (I) and the number of jets (j) in the event. We investigate a large number of cuts on the 
Pt of jets and leptons in combination with transverse sphericity, St, and missing energy, Pt- For 
clarity, we order objects by their pr, i.e. the hardest jet would be denoted j±, and we define m e s 
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and Ht as 

4 4 

m cS = Y^PT Hi) + $t, h t = ^PAxi) + $ T , (17) 



i=l i=l 

3 



where x% is a visible object and where the sum is over the first four hardest objects . In this 
analysis we define a signal that produces S events to be discoverable for a particular signature cut 
if S > max {5^,10}, where B is the number of SM background events. 

In the analysis we investigate a broad set of cuts to enhance the significance. The optimal cuts 
were found by varying the bounds on observables. First, a broad optimization was carried out where 
the varied observables include missing energy (100 GeV to 800 GeV in steps of 50 GeV), transverse 
sphericity of all visible objects (with a lower bound of 0.15 to 0.25 in steps of 0.05), number of jets 
(2 to 6 in integer steps), number of 6-jets (0 to 3 in integer steps), as well as the px of the hardest 
jet (10 GeV to 500 GeV in steps of 100 GeV) and second hardest jet (10 GeV to 250 GeV in steps of 
50 GeV). Further, this optimization includes varying the number of jets, the pt of the hardest jet, 
Pt of the second hardest jet, and the pt of all the jets. For the particular cuts, C17 and C18, that 
deal with opposite sign same flavor (OSSF) leptons, a Z-veto is applied, i.e. the leptons invariant 
mass is not in the 76 GeV to 105 GeV region. We note that in a preliminary scan, before the large 
optimization, a variation on cuts for the pt of leptons was also investigated. However, it was found 
to be of little use in enhancing the significance at low luminosity for the models we discuss. Given 
this we omitted cuts on lepton px from the large optimization. However, the tri-leptonic signal is an 
important signature for the discovery of supersymmetry via the off-shell decay of the W as well as 
other off-shell processes 861 ] . The tri-leptonic signal has been considered in the analyses at 7 TeV in 
recent works 10, lU( . A third optimization search was done in this channel by varying transverse 
sphericity (either no cut or > 0.2), missing energy (> 100 GeV, 150 GeV, 200 GeV, 250 GeV), 
number of jets (> 2, 3, 4, 5, 6), p T (ji) (no cut or > 60 GeV, 100 GeV, 150 GeV) and p T {h) (no cut 
or > 20 GeV, 30 GeV, 40 GeV, 60 GeV, 80 GeV, 100 GeV). 

A subset of cuts found using the procedure above are listed below. In choosing these cuts 
we have taken into account the uncertainty of how well missing energy can be determined in the 
early runs. For this reason we have taken lower values of $t, sometimes as low as 100 GeV. 
Better optimization can occur with other choices, specifically for larger values of pt- However, 
this requires a greater degree of confidence on how well j£y is determined in the early runs. 



3 One could also define a b jet effective mass to be the sum of the prs of the four hardest 6-jets. However, accuracy 
with which pr of the 6 jets can be determined may not be high in early runs and thus the use of &-jet effective 
mass may be experimentally challenging. Several alternate definitions of Ht appear in the literature. 
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CI: $ T > 100 GeV , 

C2: S T > 0.2, $ T > 100 GeV , 

C3: S T > 0.2, $ T > 100 GeV, n{£) = 0, PT (ji) > 150 GeV, pr{j2,h,k) > 40 GeV , 

C4: S T > 0.2, $ T > 250 GeV, n{£) = 0, p T (ji) > 250 GeV, Pt&JsJa) > 40 GeV , 

C5: S T > 0.2, $ T > 150 GeV, n{£) = 0, p T (ji) > 150 GeV, p T {j2,h,h) > 40 GeV , 

C6: S T > 0.2, $ T > 250 GeV, n{£) = 0, p T (ji) > 100 GeV, pr(j 2 ) > 40 GeV , 

C7: Sr > 0.2, ^ T > 200 GeV, n(£) = 0, p T (ji) > 30 GeV , 

C8: S T > 0.2, $ T > 200 GeV, n(j) > 2, n(^) > 2 , 

C9: S T > 0.2, $ T > 200 GeV, n(6-jets) = 1 , 

CIO: S T > 0.2, $ T > 100 GeV, n(£) = 0, n(6-jets) > 1 , 

Cll: S T > 0.2, $ T > 100 GeV, n{£) = 0, n(6-jets) > 2 , 

C12: St > 0.2, $ T > 100 GeV, n(j) > 4 , 

C13: S T > 0.2, ^ T > 100 GeV, n{j) > 4, p T {h) > 100 GeV, m cS > 400 GeV , 

C14: S T > 0.2, ^ T > 100 GeV, n(j) > 4, p T {ji) > 100 GeV, m cS > 550 GeV , 

C15: S T > 0.2, ^ T > 100 GeV, n(j) + n{£) > 4, p T (ji) > 100 GeV, H T > 400 GeV , 

C16: S T > 0.2, $ T > 100 GeV, n(j) + n(£) > 4, p T {h) > 100 GeV, F T > 550 GeV , 

C17: S T > 0.2, $ T > 100 GeV, Z-veto, n(£+) = 1, n{£^) = 1, p T (£ 2 ) > 20 GeV, p T {ji) > 
100 GeV, p T {j2) > 40 GeV , 4 

C18: St > 0.2, $ T > 100 GeV, Z-veto, n(£+) = 1, = 1, p T {£ 2 ) > 20 GeV, p T (j 2 ) > 40 GeV , 

C19: S T > 0.2, $ T > 100 GeV, ^ T > 0.2m cS , n{j) > 4 , p T (ji) > 100 GeV, 

C20: $ T > 100 GeV, n(£) = 3, p T (ii) > 150 GeV, n(j) > 2 , 

C21: $ T > 150 GeV, n{£) = 3, p T (i 2 ) > 40 GeV . 

4 In the specification of the cuts C17 and C18, the subscripts a and b indicate that they may be different flavors, 
but a Z-veto is applied only to OSSF. 
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LHC significance over channels for 1 fb of integrated luminosity 
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TABLE IX: A display of the signal significance S/ y/~B in each discovery channel for the models in Table(J| for 
1 fb -1 of integrated luminosity at the LHC. As already mentioned in Sec. (IV) for a signal to be discoverable 
we require S > max < bV~B, 10 >. 
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TABLE X: The values in the table are the integrated luminosity in units of fb -1 when the model is first 
discoverable in that channel. The table shows that many of the low mass gluino models will become visible 
with an integrated luminosity of 0.5 fb -1 , and all models become visible with an integrated luminosity of 
2 fb -1 except the model LG9 which requires an integrated luminosity of 5 fb -1 to be discovered under the 
criterion given in Sec. (IV) and in the caption of Table(IX). 
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3. Signature Analysis: Our analysis is carried out to determine the potential for discovery of 
the dark matter motivated models LG1, . . . ,LG17, for 0.5 fb _1 , 1 fb _1 , 2 fb _1 and 5 fb _1 (with a 
focus on 1 fb -1 ) of integrated luminosity with 7 TeV center of mass energy. These models exhibit 
generic features of a very broad class of SUSY models. 

We now discuss Fig.([T]), shown in Sec. (I), a bit more generally. The potential for discovering the 
models of Table (JXJ) at 1 fb _1 of integrated luminosity is given in Table (|IXp . A subset of the results 
in Table pxp are given in Fig.(fT]). Thus, Table (|IX|) displays each model's significance, S/yB, for 
the 21 cuts listed above. One finds that a particular model, LG&, often has several signatures 
that lead to large excesses of signal over background, i.e. S/VB > 5, and the set of signatures 
in which the model becomes visible varies significantly from one model to the next. Additionally, 
in studying the tri-leptonic channels, C20 and C21, we see that consistently only 6 models among 
those listed in Table dH) are discoverable (LG1, LG2, LG10, LG11, LG14, LG15) and the majority 
of the other models show less then 10 events in the tri-leptonic channels. 

The overall production cross section of superparticles is determined mainly by the production of 
squarks and gluinos. For the gluino production modes at the LHC, the cross section is determined 
by the gluino mass. However, for models with low mass gluinos, like those that appear in Tabled, 
the detectable signals at the LHC are strongly influenced by the other low- lying superparticles, i.e., 
the superparticles that are lighter than the gluino. In Fig. ([3]), an analysis of the jet multiplicity 
and the transverse momentum of the leading jets is given for LG3, LG4, LG11, LG13, and LG14. 
The distributions for jet multiplicity and jet momentum look quite different from model to model. 
For instance, the model LG11 has a gluino mass of 433 GeV, and several of its superparticles are 
lighter than the gluino, including the lighter stop and gauginos. This leads to lengthy cascade 
decay chains which produce multiple jets. Further, the mass differences between the superparticles 
in LG11 are relatively large which give rise to large momentum of the SM final states including 
jets. In contrast, models LG3 and LG4 tend to produce events with less jet multiplicity and smaller 
transverse momentum, which is due to the fact that these models having a gluino as the NLSP; 
i.e., these models are GNLSPs. For these GNLSP models the masses of the gluino and the LSP 
are correlated in the gluino coannihilation mechanism such that the mass gap is relatively small. 

Specifically, LG3 has M g — M^o ~ 50 GeV and the gluino production is, overwhelmingly, the 
dominant production mode. For this model, the gluinos decay directly to the LSP + 2 jets, i.e., 
Br(g — > {bt>Xi,qqXi)) ~ (20,80)% where q stands for first two generation quarks. We note in 
passing that generally one needs to take into account the radiative decay of the gluino, g — > gx\- 
Such a case occurs, for example, in LG4, and the decay g —> gx± dominates the branching ratio. 
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FIG. 3: (Color online) Top Left: Distribution of the number of jets without cuts. Top Right: Distribution of 
the pt of the hardest jet also without cuts. Bottom Left: Distribution of the number of SUSY events (plus 
SM background) vs. the number of jets after a cut of $x ^ 200 GeV. Bottom Right: Distribution of the 
number of SUSY events (plus SM background) vs. the px of the hardest jet after a cut of $ T > 200 GeV. 



The relatively small mass splitting in model LG3 between the gluino and the LSP (as well as the 
extreme case of LG4) makes this model class harder to discover due to the softer jets and low jet 
multiplicity, compared to other models. (For recent work on relatively small gluino-LSP splittings 
52] and Q, Q-) This feature is illustrated further in TabledTX]) and TabledH). We 
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note that in such cases where the mass gap between gluino and the LSP is extremely small, the 
effects of the ISR can be substantial for the collider signatures. We also note that although it 
can be challenging to discover events from gluino production for the GNLSP models, (depending 
on the degree of the mass degeneracy), one should keep in mind that some other subdominant 
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FIG. 4: SUSY plus SM background events vs pr(ji) at 1 fb 1 for the signature cut pr (ji) > 100 GeV, 
Pt (j'2) > 40 GeV, $ T > 250 GeV, S T > 0.2 and n(t) = for LG2, LG3, LG13. The figure illustrates the 
softness of the jets in model LG3, a GNLSP model, relative to the models LG2 and LG13. 

SUSY production modes could be detectable and become the leading signals for such models. For 
example, in model LG3, the stop is relatively light and decays entirely into a chargino and bottom 
quark, i.e., Br{t\ — > X\b) ~ 100%, and the chargino subsequently decays entirely into a neutralino 
and a W boson, i.e., Br(xt — > XiW + ) ~ 100%. Hence, this GNLSP model class may be able 
to produce discoverable leptonic events through these decay chains with upgraded center of mass 
energy and luminosity. 

Further, these features of the GNLSP models can explicitly be seen by studying the top panels 
of Fig. ([3]) and by observing the relative broadness (or width) of the pt distribution of models LG11 
and LG14 relative to LG3 and LG4. In addition, LG13, a stop NLSP model, is also peaked at low 
jet pt much like LG3 and LG4. The stop mass for model LG13 is near 200 GeV and the stop-LSP 
mass splitting is small (< 30 GeV). Thus, this model produces stops at a large rate, which decay 
(via an off-shell loop-induced and FCNC decay) into a charm quark and LSP (ii — > ex?) with a 
~ 75% branching ratio. However, the softness of jets in LG13 mimics the softness of jets in LG4 
in part due to the phase space, which explains the peaking of the distributions at low pt- The 
restricted phase space from the small mass splittings is also why the effective mass distribution for 
stop NLSPs is narrow (see j^]) relative to other cases. 

In Fig.@ we highlight the GNLSP model LG3, which satisfies the double-sided relic density 
band, along with the light stau and stop models, LG2 and LG13, which also satisfy the WMAP 
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bound via scalar coannihilations. Thus, Fig.(U]) shows jet pt, signal plus background, for the 
models LG2, LG3, and LG13 compared to the SM background alone. As discussed earlier the 
model LG3 arises from gluino coannihilations and has a relatively small mass splitting between 
the gluino and the LSP neutralino. This is to be contrasted with the model LG2, which satisfies 
the WMAP relic density band via stau coannihilations, or the model LG13, which also satisfies 
the WMAP relic density band via stop coannihilations. Because of the compressed spectra of LG2 
and LG13, there are more jets arising from the combinations of both low mass gluino and the low 
mass squark production relative to the dominant gluino production found in the GNLSP model 
LG3. This effect is exhibited in the figure. For model LG2, as the scalars are quite light, and 
even lighter than the 715 GeV gluino, the cross section for the production of squarks as well as 
the mixed squark gluino production cross sections are about an order of magnitude larger than 
the gg production. Here the gluino two body decay modes are spread out rather uniformly with 
no dominant channel. Instead the first two generation squark decay modes are short with large 
branchings. In particular, one has for the first two generation squarks, Br(qR — > Xil) ~ 100% and 
Br(qi — > X2I) ~ 32% as well as Br((qd L , q UL ) — > Xi ' + \Qu,Qd)) ~ (60 — 65)% for each decay. Thus, 
the two body decays of the first two generation squarks provide the large signal in model LG2 
even though the gluino is quite light. In addition, for LG2, the direct production of chargino pairs 
as well as chargino and neutralino production is competitive with the squark production, leading 
to leptonic decays and large lepton multiplicities. The discovery potential of the models is also 
exhibited in Fig. ([]]), where one finds that the significance of LG3 and LG4 is much less than for 
LG11. 

In Fig. ([5]) (left-panel) we show the potential for early discovery for the PAMELA compliant 
model class discussed in Table (|VIip at yfs = 7 TeV and at an integrated luminosity of 1 fb _1 . The 
displayed models have a rather large gluino production. Since the gluino, the light chargino, and 
the second heaviest neutralino are the lightest SUSY particles beyond the LSP, and the squarks are 
rather heavy for these models, the sparticle production at the LHC will be dominated by gg, Xzict 
and X1X1 production. For example, models (LG10, LG16, and LG17) have a total SUSY cross 
section of ~ (12, 15, 5) pb at leading order and the gluino production is at the level of ~ (9, 14, 4) pb, 
respectively. The chargino neutralino production makes up most of the remaining part of the cross 
section. The leading decays of the gluino are g — > xf + QQ r and g — > X2/X1 + QQ- These decays are 
subsequently followed by X2 ~ > Xi + // anc ^ Xi — * Xi + //' where /, /' are the standard model 
quarks and leptons. In particular, the lightness of the gluino in the three models (LG10, LG16,and 
LG17) gives rise to multi-jets which produce a strong signal over the background. Hence, these 
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FIG. 5: Left: SUSY plus SM background events vs m e g at 1 fb 1 of integrated luminosity for the signature 
cut p T (ji) > 150 GeV, p T (j 2j ja, k) > 40 GeV, $ T > 150 GeV, S T > 0.2 and n(£) = for the PAMELA 
compliant models. As discussed in the text LG10 is a Higgsino LSP model and LG16 and LG17 are models 
with a mixed- wino LSP. Right: The same as the left panel except for a subset of the GNNLSP models 
(with chargino and neutralino degenerate), i.e., LG7, LG8, along with the compressed models LG11, LG12, 
which in addition to a low mass gluino, also have a light stau and a light stop and have a compressed mass 
spectrum for the first two generation squarks and sleptons. 



models are good candidates for early discovery. Further, if a model of the type LG10, LG16, 
or LG17 is verified at the LHC, it would also provide a consistent explanation of the PAMELA 
anomaly. However, to fully demonstrate the validity of the models, additional luminosity would be 
needed to extract information about the neutralino mass. We do not give a detailed methodology 
for accomplishing this, but as argued in [13] it may be possible to extract information about the 
neutralino and the chargino states in the gluino decay products. 

In Fig. ([6]), we show a comparison of di-jet invariant mass distributions for the GNNLSP models 
compared to models where the gluino is positioned higher in the mass hierarchy. One sees the 
GNNLSP models (LG7, LG8, and LG10) have a relatively narrow di-jet invariant mass which 
arises from these models being dominated by the three-body decays resulting from X2X1 j X\X\ 1 
and gg production. Further, the distributions for the models LG8 and LG10 become depleted (or 
more narrow) relative to the light stau and light stop models from the three-body decays which 
result in softer jets. These subsequent decays produce an increase in the multi-jet signal compared 
to the SM background. However, the light stau and light stop models LG11, LG14, and LG15 have 
a relatively broader distribution, which arises from the compression of their sparticle spectrum. 
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di-Jet Invariant Mass (GeV) di-Jet Invariant Mass (GeV) 

FIG. 6: Left: SUSY plus SM background events vs the di-jet invariant mass (rrijj) at 1 fb" 1 of integrated 
luminosity for signature cut $ T > 200 GeV, S T > 0.2 and n(j) > 2 for the models LG7, LG8, LG10. 
Right: Same as the left plot except that the analysis is for models LG11, LG14, LG15. The left panel 
shows the light gluino models which are effectively GNNLSP models, while the right panel shows the models 
with a compressed mass spectrum for the scalars and for the light gluinos. As such the right panel shows 
distributions which are significantly broader from the squark production and decays. 

For these models, all the sparticle masses are below 700 GeV. Further, LG14 and LG15, where 
the mass spectra are compressed, the gluino is in the 31 st position of the mass hierarchy. The 
compressed spectra causes a large sampling of sparticle production, which results in a production 
of many jets with a more diverse range of momentum. 

The right panel of Fig. ([7]) shows the number of SUSY signature events plus the SM background 
in 40 GeV energy bins at 1 fb _1 of integrated luminosity vs the di-jet invariant mass for models 
LG1, LG11, and LG14. As exhibited in this figure, these models have a significantly larger di-jet 
invariant mass compared to the SM. As discussed earlier LG11 has a relatively large sparticle mass 
splittings in the scalar sector relative to the LSP mass as well as lengthy cascade decay chains 
that produce multiple final state jets with large momentum. Further, the right panel of Fig.© 
helps illustrate the effectiveness of the m e g cut. Comparing the values of C13 to C14 in Table t) LXp . 
one sees that the significance for models LG1, LG11, and LG14 increases as m e g increases. For 
the case when m eS > 400 GeV (C13) we get S/y/B = (17,34,15) and when m cS > 550 GeV 
(C14) we get S/y/B = (20,38,20), respectively, for (LG1 ,LG11, LG14). However, models LG3, 
LG10, and LG16 have a reverse effect, i.e., S/y/B = (3,23,27) for m eS > 400 GeV (C13) and 
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OSSF di-Lepton Invariant Mass (GeV) di-Jet Invariant Mass (GeV) 



FIG. 7: Left: SUSY plus background events for models LG11, LG14, LG15 vs the OSSF di-lepton invariant 
mass (m e +f-) at 1 fb" 1 for signature cut $ T > 200 GeV, St > 0.2 and n(j) > 2 with 2 leptons of any sign 
and flavor. Right: SUSY plus background events for models LG1, LG11, LG14 vs the di-jet invariant mass 
(rrijj) at 1 fb _1 of integrated luminosity for signature cut $ T > 100 GeV, St > 0.2, Pt(Ji) > 100 GeV, 
m e fi > 550 GeV and n(j) > 4. Here the peak in the distribution is a consequence of the m e g cut. 

S/y/B = (2, 15, 15) for m cS > 550 GeV (C14) for (LG3, LG10, LG16), respectively. These effects 
arise since models LG3, LG10, and LG16 have lower jet multiplicity, less missing energy, and fewer 
cascades than the models shown in the right panel of Fig. ([7]). For instance, the model LG16 cross 
section is dominated by gg production with the g dominantly decaying into x? or X\ ■ This results 
in low jet multiplicity and lower missing energy compared to models LG1, LG11, and LG14. 

4- Mass Reconstruction: We now discuss the potential to do mass reconstruction for some 
of the models with the early data. In the left panel of Fig. ([7]) we display the number of SUSY 
signature events plus background events in 20 GeV energy bins at 1 fb _1 of integrated luminosity 
vs the OSSF di-lepton invariant mass for models LG11, LG14, and LG15. The plot also displays 
the cuts used as well as the standard model background alone for comparison. In large portions 
of the figure, the SUSY signals plus the background distribution stands significantly above the 
background. The leptonic events are mostly produced from the gaugino cascade decays through 
low- lying sleptons. If the OSSF di- leptons arise from the same decay chain — > £ ± ^ =F — > Xi^^i 
the invariant mass from the reconstruction of the di-leptons obeys the following mass relations for 
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on-shell sleptons: 



m 



(18) 



In particular, for model LG15 the three sleptons, (ti, Sr, with the latter two being degenerate, 
contribute to the OSSF di-lepton events with the di-lepton invariant mass lying between Xi and 
X2- Using the sparticle masses from LG15, i.e., (M^o , , M-o ) = (121,137,240) GeV, one can 
use Eq. (fT8|) to obtain M f + f - < 92 GeV from the &rIJir decay modes. The %\ decays into Xi 
via Br(x2 &r^) — 28%, and then the right-handed slepton decays entirely into a lepton and 
Xi, i.e., Br(£^ — > x?^ ± ) — 100%. The decay of the light stau follows similarly through Br(x2 ~^ 
t^t t ) ~ 33%, and then the stau decays entirely into Xi an d a r, i.e., Br(f 1 — > Xi T± ) — 100%. 
The tau produced from the X2 decay has a subsequent leptonic tau decay with branching ratio 
of Br{r — > £u£V T ) ~ 35%. We do not attempt to reconstruct taus here, and we note that there 
are also further mixings arising from chargino decays which require flavor subtraction and other 
techniques to isolate lepton pairs coming from the same cascade decay. Further, due to the low 
statistics at the early runs of the LHC data, we do not perform a more detailed mass reconstruction 
in our analysis here. As discussed above, and as can be seen in Fig. (|7|), the mixings arising from 
other processes are rather small and the edge in the di-lepton invariant mass agrees well with the 
prediction of Eq.([T8j). 

We now discuss the reconstruction of the 6-tagged di-jet invariant mass peak. In the left panel 
of Fig. ([8]) we give an analysis of the number of SUSY event vs the 6-tagged di-jet invariant mass 
(rribb) at 1 fb _1 for the models LG1, LG7, and LG8 for the cuts displayed as well as a comparison 
to the SM background. One finds that the three models are distinguishable above the background. 
For these models, the majority of the 6-tagged di-jet events come from the gluino off-shell decay 
9 ~ * X2 + ^> which leads to an upper bound of the kinematic endpoint < Mg — M^o that 
is estimated to be in the range (300 — 322) GeV. In the left panel of Fig.(jHJ), one sees a hint of 
an endpoint forming in this region. However, for these models, the kinematic endpoint is not 
yet discernible; more luminosity would be needed, and further, additional uncertainties arise in 
the interpretation of the invariant mass endpoint due to additional cascade processes. A similar 
analysis can be given for models LG11, LG14, and LG15 in the right panel of Fig. ([8]). The source 
of the jets for the three models differ from those of the left panel of Fig.® due to their spectra. 
Further, LG11 produces a significantly larger number of jet events compared to those for LG14 
and LG15 due to its light color particles, i.e., the gluino and the stop, dominantly decaying to b jet 
final states. In addition, some of the b jets in the models from Fig.® come from the light CP even 
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FIG. 8: Left: SUSY plus SM background events vs the 6-tagged di-jet invariant mass (rribb) at 1 fb -1 of 
integrated luminosity for signature cut $ T > 100 GeV, St > 0.2, pr{jx) > 100 GeV, m e g > 550 GeV and 
n (j) > 4 for the models LG1, LG7, LG8. Right: Same as the left panel except that the analysis is for models 
LG11, LG14, LG15. As discussed in the text, there is a hint of kinematical endpoints forming for some of 
the models in the di-6-jet invariant mass plots exhibited above. 



Higgs. For example, in model LG15 Br(x2 -> Xih) ~ 30% and Br(h -»■ bb) ~ 80%. Thus with 
increased statistics one may be able to partially reconstruct events coming from the Higgs decay 
in this model and other models as well. 

More generally in Table (|X| ) we summarize the result of our analysis for the full set of integrated 
luminosities 0.5 fb , 1 fb , 2 fb 1 , and 5 fb -1 . The entries in the boxes in this table indicate the 
integrated luminosity at which a model listed in the first column will become visible in a specific 
signature channel listed in the top row. Thus, the entries in Table (pO) show that a good number of 
the models in Tabled will become visible at 0.5 fb -1 of integrated luminosity, and all of the models 
given in Table (jl]) will become visible (in at least one channel) at 5 fb -1 of integrated luminosity. 
Indeed, as discussed above, one observes that the relative mass splitting and the relative position 
of the gluino within the sparticle mass hierarchies strongly influences the discovery capability of 
the LG models. Several channels in some cases are needed to establish a signal, and the variance 
amongst channels for different models is apparent. 
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V: CONCLUSION 

In this work we discussed the sparticle landscape in the context of a low mass gluino which is 
one of the prime superparticles that has the potential of being produced as well as detected at the 
early runs of the LHC. This is due to the gluino (and also squarks) being strongly interacting and 
typically having the largest production cross section in pp collisions at the LHC for sufficiently 
low gluino masses. The low mass gluino models considered arise in a variety of settings including 
mSUGRA, nonuniversal SUGRA models, and in supergravity/string models with a very weakly 
coupled U(l) n extended hidden sector. A number of specific benchmark models were analyzed and 
found to be encouraging for discovery at both the LHC and in dark matter experiments. It is found 
that the eigencontent of the LSP in such models can vary over a wide range from the LSP being 
a pure bino, or a mixed-wino (LG16 and LG17) to the LSP being heavily Higgsino dominated 
(LG10). Further, the associated sparticle spectrum is found to be widely different with the squarks 
and sleptons being as low as 200 GeV (or even less) in mass to being significantly heavy lying in 
the (2 — 3) TeV region. The models analyzed exhibit a wide array of sparticle mass hierarchies 
and signatures. It was shown that most of the models considered will be discoverable in the early 
runs at the LHC at \fs = 7 TeV with 1 fb _1 of data while all the models will become visible with 
5 fb _1 of integrated luminosity. The models considered are consistent with the stringent bounds 
on the annihilation of neutralinos into 77 and 7Z from Fermi-LAT, and further many of the 
models considered are discoverable in the on going dark matter experiments; specifically CDMS-II, 
XENON100, and EDELWEISS-2. Further, it is found that some of the low mass gluino models 
(LG10, LG16, and LG17) can explain the positron excess observed in the satellite experiments such 
as PAMELA that probe antimatter in the Galaxy. It is shown that such models also lead to rich 
jet signatures at the LHC, thus representing a class of models which can be tested on multiple fronts. 



Note added : After the publication of this paper, CMS [87j and ATLAS [88] released their 
analyses in the search for SUSY at the LHC at -y/i = 7 TeV with 35 pb _1 . Models LG2, and 
LG(11-15) have mass scales that lie close to the observable limits quoted by ATLAS j^]. 
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